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Abstract: In future wireless networks, a significant number of users will be vehicular. One 
promising solution to improve the capacity for these vehicular users is to employ moving relays 
or car base stations. The system forms cell inside the vehicle and then uses rooftop antenna 
for backhauling to overcome the vehicular penetration loss. In this paper, we develop a model 
for aggregate interference distribution generated from moving/parked cars to indoor users in 
order to study whether indoor femto-cells can coexist on the same spectrum with vehicular 
communications. Since spectrum authorization for vehicular communications is open at moment, 
we consider two spectrum sharing scenarios (i) communication from mounted antennas on the 
roof of the vehicles to the infrastructure network utilizes same spectrum with indoor femto-cells 
(ii) in-vehicle communication utilizes same spectrum with indoor femto-cells while vehicular to 
infrastructure (V2I) communication is allocated at different spectrum. Based on our findings we 
suggest that V2I and indoor femto-cells should be allocated at different spectrum. The reason 
being that mounted roof-top antennas facing the indoor cells generate unacceptable interference 
levels. On the other hand, in-vehicle communication and indoor cells can share the spectrum 
thanks to the vehicle body isolation and the lower transmit power levels that can be used inside 
the vehicle. 
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1. INTRODUCTION 

The desire for entertainment and infotainment services 
while on board and the need for better traffic safety is 
expected to increase the traffic generated by vehicles in 
the near future [1, 2]. For the time being infotainment 
services utilize the existing macro-cellular 3GPP infras¬ 
tructure. However, in dense urban cities and/or during 
traffic jams, the capacity demand can be high and the 
existing infrastructure might not be able to support a fast 
access to the public cloud for all users. Also, nowadays, 
passengers are connected to the macro-cellular network 
directly, through the metallic vehicle body that can in¬ 
troduce losses up to about 10 ~ 20 dB [3, p.57]. In case 
metal coated windows are used, the losses can be even 
higher. Providing an adequate quality of experience for the 
users on-move requires the development of new technical 
enablers. Possible enablers can include ultra-dense network 
deployments and mounted antennas on the roof of the 
vehicles acting as moving relays [4]. 
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Nowadays, cellular deployments in urban environments 
are characterized by inter-site distances as low as 200 
m [5]. Further network densification along urban streets 
(hereafter street micro cells) can be used to increase the 
amount of load that can be served. At the same time, 
passengers inside a vehicle can connect to the micro base 
station through a gateway with an antenna mounted on 
the roof of the vehicle overcoming high penetration losses. 

The authorization scheme for vehicular communication is 
not yet standardized. One option could be to use same 
spectrum with macro-cellular network, i.e. 800 MHz to 
2.6 GHz in urban areas and either partition the spectrum 
between macro and micro layer or use the full spectrum 
under a shared spectrum access regime. Same approach 
has been proposed for overlaying indoor femto-cells to the 
macro-cellular network [6]. In this paper, we study spec¬ 
trum coexistence issues between the micro and the femto 
layer. We consider two different possibilities: (i) femto-cells 
utilize same spectrum with street micro-cells i.e. antenna 
on the roof of the vehicle generates interference indoors (ii) 
femto-cells utilize same spectrum with in-car base station 
while street micro-cells are allocated at different spectrum. 
In order to assess the impact of vehicular transmissions to 
indoor users, we develop a model for aggregate interference 
and SIR distribution. 

Modeling of aggregate interference level has received a lot 
of attention in the literature [7, 8, 9, 10, 11]. Recently, 





elements from stochastic geometry have been incorporated 
into interference related studies due to their analytical 
tractability: When the locations of transmitters are mod¬ 
eled using a Poisson point process (PPP), the Laplace 
transform (LT) of the interference distribution can be 
expressed in an integral form and for pathloss exponent 
equal to four, the integral of the inverse LT can be solved 
in a closed-form [12]. Usually, a simple power law model is 
used to describe the distance-based attenuation from each 
individual transmitter and the impact of fast/slow fading 
can also be incorporated into the model. While power law 
model is sufficient to describe distance-based propagation 
path loss in outdoor macro-cells, there are more accurate 
models to describe attenuation along street micro-cells. 
According to these models the power law changes beyond a 
certain distance breakpoint [13]. This makes the analytical 
treatment of interference distribution more involved. 

In this paper we compute the LT and the moments of 
aggregate interference generated from parked/moving cars 
at the worst case located indoor user facing a crossroad. 
We use the method of moments to approximate the distri¬ 
bution of aggregate interference. We show that the inverse 
Gamma distribution approximates closely the interference 
distribution in the upper tail. Since the outage probability 
at the indoor user is determined from the lower tail of 
the SIR distribution, modeling the interference by inverse 
Gamma distribution turns out to be a reasonable choice. 

With a low-complex and accurate interference model at 
hand we turn our attention to the SIR distribution mod¬ 
eling. We assume Nakagami-m channel model for radio 
propagation indoors and express the outage probability at 
the femto-cell in terms of the Meijer G function which in 
our cases of interests reduces to generalized hypergeomet¬ 
ric function widely available in most of today’s software 
packages for numerical computations. By using it, we 
study how the density of cars, the uplink transmit power 
level and the car isolation impact the outage probability at 
the femto-cell. We show that mounted roof-top antennas 
generate unacceptably high interference level at the femto- 
cell while in-car and femto-cell communication can coexist 
in same spectrum. Also, we utilize the proposed model to 
illustrate that the most significant part of the interference 
comes from the street facing the femto-cell. Based on this 
observation, we argue that for mounted roof-top antennas 
frequency planning between street micro-cells and indoor 
femto-cells could be an option for performance enhance¬ 
ment. 

2. SYSTEM MODEL 

We consider a dense urban city with street layout following 
the Manhattan grid. The building block length is denoted 
by d and the street width is denoted by r. Let us define the 
distance between neighboring streets as D = d\r , see Fig. 1. 
Along a street, the vehicles are distributed according to a 
PPP model. At particular snapshot, each vehicle is active 
in the uplink with a certain probability. As a result, the 
distribution of uplink transmissions along a street follows 
a PPP too. We denote the PPP on the fc-th vertical street 
by with density (i.e. number of uplink transmissions 
per meter) A& and on the k! th horizontal street by 
with density A k'■ 



Fig. 1. System illustration. 

The uplink transmit power level is denoted by P t and the 
fading channel h for vehicular transmissions is assumed to 
be Rayleigh with mean equal to unity, E{h} = 1. The 
distance-based pathloss in Manhattan grid, l(x 1 ,^ 2 ); is 
a function of the NLOS distance, x±, as well as of the 
LOS distance x%, see also Fig. 1. For a vehicle located 
at coordinates (xi,X 2 ), the distance-based pathloss to the 
origin can be read as Z(xi,X 2 ) = Cx]~ 2 x^“ where the 
attenuation exponent a along the LOS is taken equal to 2 
and typical values for the NLOS are 4^6. One may notice 
that for the k -th vertical street X 2 t k = kD, k > 0 while for 
the vertical street facing the considered femto-cell, k = 0 , 
we have X 27 - = 1- Also, the parameter C accounts for 
antenna gain, car isolation (in the case of in-car base 
station) and wall losses with typical wall attenuation in 
the range of 15 ~ 20 dB. 

We want to identify whether uplink vehicular commu¬ 
nication and indoor femto-cells can coexist in the same 
spectrum. In order to do that we evaluate the outage 
probability for given SIR target 7 at the worst-case lo¬ 
cated femto-cell facing a cross-road, see Fig. 1. In order 
to evaluate the SIR distribution we need a model for 
the aggregate interference distribution and a model for 
the wanted signal power level at the femto-cell. Assuming 
Nakagami-m fading for indoor propagation, the wanted 
signal power is Gamma distributed with shape m and scale 
P rx /m where P rx denotes the mean received signal power. 
Next, we propose a model for the interference distribution 
and evaluate its moments. 

3. MODELING INTERFERENCE DISTRIBUTION 

The interference level can be expressed as the sum of 
interference levels from all active car base stations on all 
vertical and horizontal streets. 

k k' 

=ypM(x hi ,x 2 ,k)+y, Pthyiix^ji,x 2 ,k>) 

k k' 

where xi 7 and X 2 ,k are vertical NLOS distance and hori¬ 
zontal LOS distance, respectively, between *-th car on fc-th 
vertical street to the crossing toward the victim building, 
and x\y and X 2 ,k' are horizontal NLOS distance and 
vertical LOS distance, respectively, between i'-th car on 
fc'-th horizontal street to the crossing toward the victim 
building. 













































Based on the assumption that vertical/horizontal streets 
in Manhattan two-dimensional grid are symmetric and a 
constant car density is used on all the streets, we have 
same density for the processes 4>fc/Vfc, k', or, A& = 
A k 1 = Wk,k r . Also, the independent property of PPP 
allows us to focus only on one type of street (hereafter, 
vertical streets) and incorporate the impact of horizontal 
streets using a scaling factor equal to two. Next, we 
characterize the distribution of aggregate interference from 
uplink transmissions located on the £;-th vertical street, 
I $ k , in terms of its LT. 

Lemma 1. For Rayleigh fading channels, i.e. hi ~ exp(l), 
the LT of the aggregate interference from car base stations 
distributed as one-dimensional PPP on the fc-th vertical 
street in a Manhattan grid is 


£j **«( s ) = exp 



( 1 ) 


where 2 = P t C and /3k = x% k . The proof is given in 
Appendix. 

Proposition 2. The total interference from all vertical 
streets in a Manhattan grid is 


P'1, 


*00(s) = ex P _7rA { 1 + < ^(f) jD *}Vsz (2) 


where £(•) denotes the Riemann zeta function. The proof 
is given in Appendix. 


Note that the integral range (refer to the Proof) used in 
Eq.(l) is from 0 to infinity, i.e. singular pathloss model. 
This range causes inaccurate interference modeling in the 
near-field, since the interference level can become infinite 
when X\ and X 2 ,k become zero. In a more practical case, 
non-singular path loss model is taken into account by 
limiting X\ > 1, and apply the integral J) ( -)dxi . The 
results given in Eq. (1) and Eq. (2) can be regenerated 
through the following Lemma and Proposition. 

Lemma 3. For the non-singular path loss model, the inter¬ 
ference from the fc-th vertical street in a Manhattan street 
is 



0) = ex p 

The proof is given in Appendix. 

Proposition /. For the non-singular path loss model, the 
total interference from all vertical streets in a Manhattan 
grid is 


£/#(*> ( s ) = ex P 


v—■v / SZ 

-2A > 4 / — arctan 

k—Q 



(3) 


The integral for the inverse LT in Eq. (3) does not exist 
in closed-form. As a means to provide simple and useful 
expressions, the use of approximations for the aggregate 
interference is motivated. Since the LT of the interference 
with the singular pathloss model in Eq. (2) resembles 
the Characteristic function of the Levy distribution, we 
approximate the interference by a fitted inverse Gamma 
distribution. 


The first two moments of the interference can be found 
based on the LT E(/|,) = (—l) l ^ T £/ # (s)| g=0 and ex¬ 
pressed as following property of the 


E{U}=2\J2 3-=2Az 

fc =0 p k 


1 + 


4A 


OO 9 

z z 


e {4}= t Et?2+ e U4 2 = 


CW 

D a 
4A 2 2 


1 + 


C(2a ) 

D 2o 


-E{/ $ } 2 . 

fc=o^fc ' ■ - 

The fitted inverse Gamma distribution has scale a = 
+2 and shape b = E{/$}(a — 1). The inverse 

of the interference distribution is a Gamma distributed 
random variable with scale 
Gamma(a, 1 /b). 


and shape 1/6, /, 
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4. MODELING SIR DISTRIBUTION 


For the performance evaluation of wireless systems, assess¬ 
ing the quantity of the signal to interference plus noise 
ratio (SINR) is required. Based on our system model 
assumptions, the performance at the worst case femto-cell 
is limited by the interference level and the impact of noise 
power can be ignored. The SIR level is 

SIR=^L 

where P rx is the wanted received signal level. 

For ensuring satisfactory services to users of the femto-cell 
network, the outage probability (i.e. probability that the 
SIR value is below a target SIR) must be maintained under 
specific value. The outage probability can be computed 
using the following Lemma. 

Lemma 5. The outage probability when the source desti¬ 
nation fading is Nakagami-m distributed, with mean signal 
level P rx is 


LTL 1 / j. \ j 7 2 

P(S/R<7) = l-E i /f- 5? £,„.,«) (4) 

where £ = =p -. The proof is given in the Appendix. 

Eq. (4) indicates that the LT of the aggregate interference 
generated by the car base stations plays a central role 
in quantifying the outage probability in the femto-cell 
network. While the outage probability is completely char¬ 
acterized by the LT of the interference, it is also useful to 
approximate the SIR distribution by some known function 
so as to assess its mean and higher moments in a low- 
complex manner. 

Note that the useful signal level at the femto-user in 
the Nakagami-m fading channel is Gamma distributed 

with scale m and shape 0 = Prx ~ Gamma(m,0). 
Hence, with I^ 1 ~ Gamma(a, 1/6) presented in Section 3, 
the SIR, P rx ■ /,/ 1 , can be expressed as product of two 
independent Gamma random variables. By normalizing 
respective shape parameters, the normalized SIR j n = | ■ 
becomes the product of two Gamma random variables 
with shape equal to unity, gi = P rx S l ~ Gamma(m, 1) 
and g 2 = L $6 ~ Gamma(a, 1) and 7 n = gi ■ <72 • In that 
case the PDF of the normalized SIR distribution can be 
expressed in terms of the Meijer G function[14] 

p/ , ^02 ( 771,0 lo'ra) 2777 , 2 • Km— a(2^/7ra) 

l7nj = r(m)r(a) ~ r(m)r(a) 



















Fig. 2. Approximations of interference level distribution 
using Gamma and Inverse-Gamma distribution with 
non-singular path-loss model, when (a) the interfer¬ 
ence from the closest street k = 0 is not considered 
and (b) the interference from all vertical streets is 
considered. 

where T(-) is the gamma function, G(-) is Meijer G- 
function and AT m _ a (-) is the modified Bessel function of 
the second kind with order (m — a). 

The CDF of the normalized SIR distribution is derived 
using integral properties [14] of the Meijer G-function, 


(b \ GlKm^a, 0 \ln) 

\e'u~ ln ) r(m)r(a) 


(5) 


where G\ 3 ( m ^ 0 | 7 „) 
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i F 2 (to; to +1 — a, m + 


1 ; In) + 1 ' — a u 1 " 1 f 2 (a; a + 1 - TO, a + 1; j n ) and p F q (-) is 
the generalized hypergeometric function for integer p and 


q- 


Eq. (5) can be used as the CDF of the SIR distribu¬ 
tion by properly shifting the axis P ( | • jb < 7 „) = 


P (if<!l")= P (A<l)- 


5. NUMERICALS 


Fig. 3. Comparison of outage probability for simulation 
result and analytical result in Eq. (4) with respect to 
car base station density A and car isolation p when 
Nakagami fading parameter to = 1, mean wanted 
signal level is P rx = —40 dBm and SIR target is 
7 = 15 dB. 



Fig. 4. SIR distribution when car base station density 
is A = 0.1 Cars/m, the mean wanted signal level is 
Prx = —40 dBm, Car isolations are (a) r) = 0 dB and 
(b) rj= 20 dB 


In order to verify the suitability of interference and SIR 
models presented in Section 3 and Section 4 we consider 
cars deployed in a city centre which is modeled as a square 
with side equal to 2 km. Street width is taken equal 
to r = 20 m and building block is a square with side 
d = 50 m. For NLOS propagation the pathloss exponent 
is taken equal to a = 4. The transmit power level in the 
uplink is P t = 100 mW, the wall attenuation is 15 dB, 
the attenuation constant is 10” 3 . For mounted roof-top 
antenna, car isolation is equal to r/ = 0 dB and for in- 
car base station rj = {10,20,30} dB are considered. For 
indoor propagation we present results for Rayleigh fading, 
to = 1. The mean wanted signal level at the femto-cell 
is taken equal to P rx = —40 dBm. In order to evaluate 
the outage probability at the femto-cell we assume a SIR 
target equal to 7 = 15 dB. 

Firstly, we check the accuracy of fitted inverse Gamma 
distribution for approximating the aggregate interference. 


For the parameter settings presented above the scale and 
the shape of the fitted distribution are a = 2.3 and 
b = 8.2 • 10 -9 respectively. Recall that the inverse Gamma 
approximation was selected on the basis of the LT of 
the interference for singular pathloss model in Eq. (2) 
which resembles the characteristic function of the Levy 
distribution. For comparison purposes we illustrate the 
accuracy for a fitted Gamma distribution too, see Fig. 2. 
Both distributions are accurate in the upper tail which 
would determine the accuracy of the approximation in the 
lower tail of the SIR distribution. However, the inverse 
Gamma distribution achieves better approximation over 
the full distribution body. In Fig. 2, it is interesting 
to note that the mean interference level decreases by 
approximately 50 dB if the impact of the street facing the 
femto-cell of interest, k = 0 , is not considered. 

Secondly, we utilize Eq. (4) for approximating the outage 
probability at the femto-cell as a function of the density 































A of uplink transmissions, see Fig. 3. With car isolation 
equal to 20 dB or higher, the outage probability decreases 
at acceptable values even for a high density of vehicles. 

The outage probability has been approximated by evalu¬ 
ating the LT of the interference distribution and also its 
derivatives at certain point, see Eq. (4). The amount of 
computations might be high particularly when the shape 
m of the Nakagami distribution is high. The CDF of the 
SIR distribution has been approximated in Section 4 using 
the Meijer G function. The approximating CDF can be 
useful for estimating not only the outage probability but 
also the moments of the SIR distribution. In Fig. 4 one 
can observe that the Meijer G function is accurate in the 
lower tail. However, the mean SIR is underestimated by 
2.5 dB. 

6 . CONCLUSION 


pathloss model l(: r, X2) = Cx 2 x 2 “ and (p 6 ) follows from 
the integration rule / Q °° -^p-^dx = [15, 2.172], 

Appendix B. PROOF OF PROPOSITION 2 

Since the PPP:s along different vertical streets are in¬ 
dependent among each other, the LT of the aggregate 
interference from all vertical streets is equal to the product 
of the LT:s from each vertical street 

k 

(?i) e -wA(n-C(f )D~Tt )Vii (gj) 

where (pi) follows from the Riemann zeta function, £(s) = 
EZik~ s [15, 0.233], 


In dense urban cities vehicular communication may use 
same spectrum with indoor femto-cells. As a result, better 
quality of experience for the end-users while on board 
would most probably come at the cost of higher inter¬ 
ference generated indoors. We developed a model that is 
useful for assessing the outage probability at the femto-cell 
due to the interference generated from car base stations. 
With mounted antennas on the top of the vehicles the 
outage probability becomes prohibitively high given that 
the density of vehicular transmissions is high too. On the 
other hand, with in-car base stations, the isolation due to 
the car shell and the possibility to use lower power levels 
inside the car make it possible to maintain a low outage 
probability at the femto-cell even for a high car density. 
However, in that case, additional spectrum resources are 
needed as only in-car communication uses same spectrum 
with indoor femto-cells while street micro-cells are allo¬ 
cated at different spectrum. 

APPENDIX 


Appendix A. PROOF OF LEMMA 1 

The distribution of the random variable I$ k is character¬ 
ized in terms of its LT which is given by 


= E[e—'**] =E 
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Appendix C. PROOF OF LEMMA 3 


For the non-singular pathloss model, the interference from 
fc-th vertical street is derived by following similar steps to 
Lemma 2. Only the lower integration limit is different. 
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where (pi) follows from the integration rule da; = 

yf-arctan(^) [15, 2.172], 


Appendix D. PROOF OF LEMMA 5 


P(SIR > 7 ) = P (hP rx //$ > 7 ) = P (h > 7 I<s>/P r x) 
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where x = =3— and £ = mx. Equality (pi) follows from 

denoting // # (t) the PDF of the interference distribution, 
(p 2 ) is based on the fact that when the source-destination 
fading is Nakagami-m distributed, then the CCDF of h 

is P (h > x) = e~ mx , and (p3) holds true 

according to the following Laplace transform property 

/ 0 °° e~ st t i fi^{t)dt= (-iy-£rC{fii.}(s). 


where (pi) follows from the i.i.d. distribution of the fading 
h, (p2) follows from the generating functional [7] of the 
one-dimensional PPP, (p3) follows from the symmetry 
of the interference around the origin, and (p4) follows 
from the fact that the fading h follows an exponential 
distribution with mean equal to unity. (p5) follows the 
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